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■ Effects of CP violation on the supersymmetric electroweak correction to the 

oo ; 

p parameter are investigated. To avoid the EDM constraints, we require that 
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are 0) < ICT 2 and the non-universal trilinear couplings Af = (0, 0,Aq) and 
also assume that gluinos are heavier than 400 GeV. The CP phase (fit = 



arg(^4o) leads to large enhancement of the relative mass splittings between I2 
and bi(ti), which in turn reduces the one- loop contribution of the stop and 

X ' 

j_j ■ sbottom to Ap. For small tan/3, such a CP violating effect is prominent. We 

also study how much the two-loop gluon and gluino contributions are affected 
by the CP phase. Possible contributions to the p parameter arising from the 
Higgs sector with CP violation are discussed. 
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The minimal supersymmetric standard model (MSSM) is the best motivated extensions 
of the standard model (SM). As is well known, in the MSSM, there are many new sources 
of CP violation beyond the Cabibbo-Kobayashi-Maskawa (CKM) phase, arizing from the 
complex soft supersymmetry (SUSY) breaking terms, i.e., the Majorana gaugino masses Mi, 
the trilinear scalar A terms, the bilinear scalar B terms, as well as from the /x parameter 
which is the bilinear mixing of the two Higgs chiral superfields in the superpotential. While 
those supersymmetric CP violating (CPV) phases could give large contributions to the 
neutron/electron electric dipole moments (EDMs), they are strongly constrained by the 
current experimental measurements of the neutron/electron EDMs, which for the neutron 
is d n < 1.1 x 10~ 25 e cm and for the electron is d e < 4.3 x lCT 27 e cm ||. To resolve 
this problem, several possibilities have been suggested. The first one is to make all phases 
small (O(10~ 2 )) [^]. But, such small phases would require a significant amount of fine-tuning 
and are thus undesirable. The second is to take the SUSY spectrum heavy in the several 
TeV range which lies outside the reach of the accelerators Another option suggested 
recently is to use the cancellations among the various contributions to the neutron/electron 
EDMs, allowing for large CPV phases and a supersymmetric spectrum that is still within 
the reach of the accelerators ||. Finally, an interesting possibility is to take a slightly 
non- universal scenario for the trilinear couplings Af ||[7|]. As shown in Ref. |7], requiring 
that arg(/i) < 10~ 2 and Af = (0,0, Aq), and assuming that gluinos are heavier than about 
400 GeV ||, one may significantly reduce the size of the neutron/electron EDMs due to 
Weinberg's three-gluon operator || well below the current experimental limit. 

Another CP violation associated with the Higgs boson sector in the MSSM can come from 
the one-loop corrections of the MSSM Higgs potential through soft CPV Yukawa interactions 
involving squarks, i.e., A terms. As shown in Ref. an immediate consequence of Higgs 
sector CP violation in the MSSM is the generation of mixing mass terms between the CP 
even and CP odd Higgs fields. Thanks to those scalar-pseudoscalar mixing mass terms, 
the small tree-level mass difference between the heaviest Higgs boson and the CP-odd scalar 
may be lifted considerably. The phenomenological consequences of Higgs sector CP violation 
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have been studied P,|TO||. 

A possible way to probe SUSY is to search for the virtual effects of SUSY particles which 
would be sizable enough to be detected in the present experiments. In most experimentally 
accessible processes, the heavy SUSY particles decouple from the low energy electroweak 
observables. However, similar to the SM, if there happens large "custodial" SU(2)y break- 
ing, the electroweak observables which are parameterized by the p parameter defined as 
the ratio of the strengths of neutral and charged currents at vanishing momentum transfer, 
can severely be affected. In the MSSM, the main potential source of SU(2)v breaking is 
the splitting between the stop and sbottom masses. In addition, there can be large splitting 
within a supersymmetric multiplet, which also affects the p parameter. From the numerical 
analysis, it is well known that the leading contribution of the squark loops, in particular stop 
and sbottom loops, to the p parameter can reach at the level of a few 1CT 3 which lies within 



the range of the experimental sensitivity [T^| . There are also small additional contributions 



coming from the exchange of the additional Higgs bosons, the charginos and neutralinos. 
Some discussions of supersymmetric contributions to the p parameter already exist in the 
literature []I2UI3| . Moreover, higher order corrections have also been studied so as to treat 
the SUSY loop contributions to the electroweak observables at the same level of accuracy as 
the SM contribution [Q. However, so far, the analyses have been done within the context 
of the CP conserving (CPC) MSSM. Since new CPV phases may affect the SUSY spectrum 
and generate additional couplings which do not appear in the CPC MSSM, it would be quite 
interesting to study the effects of the CPV SUSY phases on the electroweak observables via 
the p parameter. 

In this respect, the purpose of this letter is to examine how much the p parameter can 
be affected by the non-vanishing CP phases in the MSSM. In this analysis, we will impose 
the universal conditions on the soft SUSY breaking terms which provide three complex 
parameters, i.e., universal gaugino mass mi/2, the trilinear soft term Aq and the bilinear 
soft term B . In addition, there is a complex mass parameter p. However, not all phases 
of the four complex parameters are physical fll5 |. It is always possible to make a phase 



transformation on the gaugino fields so as to make m^j real. From minimization conditions 
on the Higgs potential, one can make the phase of B equal to the phase difference of the two 
Higgs doublets in the MSSM. As a result, there are only two independent CPV phases in 
the MSSM with universal soft SUSY breaking terms, which are usually chosen to be arg(Ao) 
and arg(/i). To satisfy the EDM constraints, we will simply take arg(/x) to be zero and the 
trilinear coupling of the Higgs bosons to the squarks of the first and second generation to 
be much smaller than the one of the third generation, as mentioned above, and assume that 
the gluino mass to be larger than 400 GeV. Thanks to this non- vanishing phase arg(Ao), 
as will be shown later, the stop spectrum and their mixing angle are changed, which lead 
in turn to the shift of the contributions to Ap. Also, since the neutral Higgs boson mass 
eigenstates in the CPV MSSM are mixtures of CP even and CP odd states, new gauge-Higgs 
couplings are induced at the tree level through those mixings. These couplings can generate 
the additional radiative corrections to the gauge boson self-energies, and may thus affect the 
p parameter. However, the contributions coming from the gaugino and higgsino exchanges 
will not be affected by the non-vanishing CPV phase arg(A ) in the leading order because 
they are affected by only the CPV phase arg(/x) which is neglected in this analysis. Thus, we 
shall not consider those contributions. In addition, we shall study how much the dominant 
two-loop contributions via gluon and gluino exchanges can be shifted by the effect of the 
CPV phases. The numerical calculation will be done with the help of the program package 



FEYNHIGGS |16J which is modified so as to calculate the CPV SUSY contributions to the 
p parameter. 

The p parameter is represented by 

i a _ n zz (o) iw(o) 

P l-Ap' np Mf Ml ' [L) 

where H.ww(zz)(Q) denotes the transverse parts of the W/Z-boson self-energies at zero mo- 
mentum transfer. In the SM, the main contribution comes from the large splitting between 
the top quark and the bottom quark masses and is given by 

Ap 5M = ^^oK,mg), (2) 



where N c is the color factor and the function Fq is given by 

, > _ x , > 

F (x,y) = x + y In-. (3) 

x-y y 

In the limit of large mass splitting it becomes proportional to the heavy quark mass squared, 
i.e., F (mf,ml) ~ m 2 . Including QCD corrections, the SM contribution of Ap is Apf M = 
-Ap$ M fe) (1 + tt 2 /3) 0. For a s (M z ) ~ 0.12, the QCD correction reduces the one-loop 
result by about 10% and shifts m t upwards by about 10 GeV |I8|. 



(A) Squark contributions: Since the scalar partners of the light quarks are expected 
to be almost degenerate in most supersymmetric theories, their contributions to Ap are 

negligible and thus only third generation will contribute. The stop mass matrix is given by 

/ mf + m 2 + D L m t m* LR \ 

M i = \ 2 ( 4 ) 

V m t m LR mf R +m z t +D R J 

where mm = A — p*/ tan/5, Dl = (| — § sin 2 0y/} cos 2/5M| and Dr = | sin 2 9w cos 2(3 M|. 

The stop fields (£l,£,r) are linear combinations of the mass eigenstates tip so that ti = 

cos 6*^1 — sin^-e - ^*^, in, = sin^e'^fi + cos%t 2 , where 

mf - mf +D L -D R 
and the mass eigenvalues of stop fields are given by 

1 



mf 

ti.a 



2 



mf + mf + 2m^ + D L + D 



n 



T - m| B + £> L - £> fi ] 2 + 4m 2 |m i/? | 

l m Li?J = |^4o ~ A 4 */ tan /3| 



1/2 



(6) 



^o| 2 + |p| 2 /tan 2 /3 - 2|A o ||p|/tan/3cos0 t , (7) 

where <pt = arg(y4o). As expected, the non- vanishing phase <pt will change the stop mass 
spectrum which can in turn lead to the shift of Ap. Neglecting the mixing in b sector, the 
contribution of (t, b) to Ap is presented at one-loop order by: 

Apfb = 8V%r 2 ^~ Sin2 9l C0S2 9iF ° ^ ' ^ 

+ cos 2 f ~F o (m£ , m} L ) + sin 2 %F (m|, mj- )], (8) 



Before examining how much the value of Ap^ can be affected by the phase (f> t , let us 
first investigate the mass splittings between the two squark mass eigenstates which may 
indicate to what extent the contribution of (t, b) to Ap is deviated. In Fig. 1 a (b), we 
present numerical estimation of the relative mass splittings between the two squark mass 
eigenstates, jm^ — m~J/m~ and |m| — ra~ |/m~, as a function of the phase <f> t , for tan/3 = 
1.6(20), m q = 200 GeV, \A \ = 200 GeV and |p| = 100 GeV. Here we take < fa < %. As 
the phase <p t increases, the both relative mass splittings increase. The effect of the maximal 
CP violation corresponding to <f> t = n may lead to the enhancement of the relative mass 
splittings by about 90%(40%) for |m| — m~J/m~ (\m~ 2 — \/m~) compared to the CPC 
case. This is attributed to the fact that the absolute value of m LR monotonically increases 
with the increasing phase <p t , which makes the mass splitting between tiipL) and I2 larger, 
and the mixing angle 61 closer to 7r/4. Note that the mixing angle 9f close to 7r/4 are naturally 
obtained because the off-diagonal elements of the stop mass matrix may be larger than the 
difference of the diagonal entries. As one can see from Fig. 1 a, in the case of small tan (3, 
the splitting |m| — m~J/m~ is larger than |m| — m~J/m~ for (f> t = (CPC case), whereas 
the latter becomes larger than the former for large <p t . Thus, the splitting originating from 
the left-right mixing of stops may be more important for the CPV case with large phase. In 
the case of large tan f3, jm^ — m-J/m- is larger than |m| — m~J/m~ for all 4>t, and the mass 
splittings are weakly dependent on <f> t since the term concerned with the phase 4>t in \ m LR\ 
is suppressed by the large tan/5. We also observed that the mass splitting \m~ — rn~ bL \/m 2 ~ 
is much smaller than the other two. For the same input values of the SUSY parameters, the 
dependence of Ap^ on the phase <p t is shown as solid lines in Fig. 2. Although the CPV case 
provides larger mass splitting than the CPC case, Ap^ is monotonically decreased with the 
increasing <f> t . The effect of the CPV phase reduces the value of Ap^ for the CPC case by 
about up to 35%. The reason is that as the phase increases, the mixing angle d~ t gets close 
to the maximal mixing and F (m|,m|) increases faster than F (m|,m^) in Eq. (8) so 
that they lead to destructive contribution to Ap^. As \ A§\ increases, Ap^ tends to decrease, 
whereas it tends to increase with the increasing We have observed that the effect of the 
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CPV phase is diminished as long as rriq becomes larger than \Aq\ and \p\. This is because 
m t 1 ~ m t 2 ~ m 6 L i n the limit of large m,. 

The two-loop QCD correction to Ap^ is given by Eq.(8) in Ref. ||14]| , and the two- loop 
contribution mediated by gluino exchange is also represented in Ref. |[L4|| . Those gluon and 
gluino contributions add up to about 30% of the one-loop value in the CPC case. 

In Fig. 3, the dependence of the dominant two-loop QCD correction and gluino contri- 
bution to the p parameter on the phase <f>t is plotted by the dashed line and dotted line, 
respectively, for tan/3 = 1.6 (a) and tan/3 = 20 (b). As the phase <pt increases, the two-loop 
gluonic SUSY contribution decreases, whereas the gluino contribution increases. The reason 
that the gluino contribution increases with the increasing phase <p t is that its contribution 
depends inversely on the mass splittings between ti(b£) and t^- In particular, for small 
tan/3, the value of the gluino contribution becomes larger than that of the gluonic SUSY 
contribution at large 4>t- 

(B) Higgs sector contributions: There are also possible contributions arising from 
the Higgs sector to the p parameter. In the CPC case, it is known [12] that the Higgs boson 
masses and couplings to gauge bosons are related in such a way that they lead to large 
cancellations in their contributions to the p parameter, which is at most of order of 10~ 4 . In 
the decoupling limit where the heavy neutral CP-even Higgs, CP-odd Higgs and the charged 
Higgs bosons are nearly degenerate and their couplings to gauge bosons tend to zero, while 
the lightest CP-even Higgs boson reaches its maximal mass value and has almost the same 
properties as the SM Higgs boson. Then, the contribution of the Higgs sector of the MSSM 
to the p parameter is practically the same as in the SM. However, while the contribution of 
the SM Higgs sector gives rise to logarithmic log M^jMz in the limit of large Higgs mass 
which may reach about 10~ 3 order, the MSSM contribution reaches at most 10 -4 due to the 
upper bound on the mass of the lightest Higgs boson. 

Now, let us take into account the contributions coming from the neutral Higgs sector 
with CP violation. The decomposition of the two Higgs doublets is given by 
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$i=l 1 I , <5> 2 = e ie \ - I, (9) 

V + 0i + «ai)/ \^(w2 + 02 + «a 2 )/ 

where i> i and u 2 are the vacuum expectation values and 9 is their relative phase. As shown 



in Ref. ||19|| , the Higgs quartic couplings receive significant radiative corrections from Af 
terms. The vacuum expectation values t>i,t> 2 and the phase 9 are determined from the 
minimization conditions on the Higgs potential of the MSSM. The CP odd fields are rotated 
to ai = cos /3G° — sin Pa and a 2 = sin (3G° + cos [3a so that the Higgs potential shows up a 
flat direction with respect to the Goldstone field G°. The no n- vanishing CPV phases in A t 
and Ab can lead to the non- vanishing 9, for which the scalar-pseudoscalar mixing mass terms 
are generated and the neutral Higgs boson mass matrix in the weak basis (G°, a, </>i, 2 ) can 



M 2 M 2 S 



be written by 

Ml = I ^ • ' (10) 

M 2 SP Ml 

where Mj, and Mj are 2x2 matrix for CPC parts in the basis (G°, a) and (4>i, </> 2 ), respec- 
tively, and M 2 PS = (M| P ) T describes the CPV mixings between (G°,a) and (0i,0 2 ). From 
the tadpole condition, the Goldstone field G° does not mix with the other neutral fields and 
thus becomes massless. Then, the neutral Higgs mass matrix Mq reduces to a 3 x 3 matrix 
M 2 , which is spanned in the basis (a, 1; 2 ). The Higgs boson mass matrix M 2 can be diago- 
nalized with the help of an orthogonal rotation matrix U: U T M 2 U = diag(Mfj 3 , Mj^, M^J. 

Then, the contributions of the Higgs sector with CP violation to the p parameter can be 
given by 

£(1 - U 2 u )F (M 2 Hi ,M 2 Hc ) -± j2(UuU 2j - U 2i Uij) 2 F (M 2 Hi ,M 2 H .) ) (11) 

,i=l A hi / 



X 



where Mh c denotes the charged Higgs boson mass and 



Ap s H M (M 2 ) = 16? J M 2 [Fo(M 2 ,M^)~F (M 2 ,Ml 



a 

+ 



IGtts 2 



w lv± w 

M 2 M 2 1 M 2 M 2 

M 2 -M 2 W ° g Ml~^M 2 -M 2 ° S Mj 



w 



(12) 



corresponds to the SM Higgs boson contribution to Ap. 

In Fig. 3, we present the dependence of the Higgs boson sector contribution to Ap on 
the phase (fit (solid lines). Similar to the one- loop contribution of stop quark to Ap, the 
one-loop contribution of Higgs sector is decreased with the increasing phase <fi t . This is 
mainly because the relative mass splittings among three neutral Higgs bosons decrease as 
the phase (fit increases as shown in Fig. 4. In the case of small tan f3, |M| — Mf |/Mjj (solid 
line) is larger than |M| — M\\jM\ c (dotted line). For only small tan/3, the contribution of 
the Higgs sector is larger than the two-loop contributions mediated by the gluon and gluino 
exchanges, and those three contributions are almost the same at (fit = it, which increase the 
one-loop contribution of stop and sbottom by about 30 — 35%. For large tan (3, the Higgs 
sector contribution to Ap is negligibly small. 

In summary, effects of CP violation on the supersymmetric electroweak correction to 
the p parameter have been investigated. To avoid the EDM constraints, we have required 
negligibly small arg(p) and the non-universal couplings Af — (0,0, A ) and also assumed 
that the mass of gluino is larger than 400 GeV. The CP phase (fi t leads to large enhancement 
of the relative mass splittings between t 2 and b L {ti), which in turn reduces the one- loop 
contribution of the stop and sbottom to the p parameter. For small tan/5, such a CP 
violating effect is prominent. We have also studied how much the two- loop gluon and gluino 
contributions are affected by the CP phase. Possible contributions to the p parameter arising 
from the Higgs sector with CP violation have been discussed. 

The work of JDK is supported in part by the Korea Institute for Advanced Study. 
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FIG. 1. The relative mass splittings between the two squark mass eigenstates, |m| — m| |/m| 
(solid line) and |m| 2 — m~ (dotted line) as a function of the phase <fi t for tan/? = 1.6 (a) and 

tan/3 = 20 (b). We take m~ q = 200 GeV, |A | = 200 GeV, and = 100 GeV. 
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FIG. 2. The contribution of stop and sbottom to Ap as a function of the phase <pt for the 
same input values of the SUSY parameters as those in Fig. 1. 
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FIG. 3. The contributions of the two-loop QCD correction (dashed line), the two-loop gluino 
(dotted line) and the Higgs boson sector (solid line) to Ap as a function of the phase <pt for 
tan/3 = 1.6 (a) and tan/3 = 20 (b). The input values of the SUSY parameters are taken to be the 
same as those in Fig. 1, and the gluino mass and the charged Higgs boson mass to be m g = 500 
GeV M Hc = 200 GeV, respectively. 
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FIG. 4. The relative mass splittings among three neutral Higgs bosons: |M| — M| (solid 
line) and |M| — M^\/M^ (dotted line). The input values of the SUSY parameters are taken to 
be the same as the case of Fig. 3. The lightest (heaviest) Higgs masses are denoted by M2{M^). 
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